Anomalous thermal expansion in a-titanium 
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We provide a complete quantitative explanation for the anisotropic thermal expansion of hep 
Ti at low temperature. The observed negative thermal expansion along the c-axis is reproduced 
theoretically by means of a parameter free theory which involves both the electron and phonon 
contributions to the free energy. The thermal expansion of titanium is calculated and found to be 
negative along the c-axis for temperatures below ~ 170 K, in good agreement with observations. 
We have identified a saddle-point Van Hove singularity near the Fermi level as the main reason for 
the anisotropic thermal expansion in a— titanium. 
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The most general aspects of the chemical bonding 
in the transition d metals can be understood from the 
Friedel model [ij, explaining the trends in equilibrium vol- 
ume, bulk modulus and cohesive energy. The transition 
metals are found to crystallize at low temperatures in 
the cubic fee and bcc structures, and the hexagonal hep 
structure , which can be qualitatively explained from a 
band filling of itinerant d-states 3]. In addition, the De- 
bye model reproduces the thermal volume expansion with 
a rather good accuracy Hence, with a seemingly good 
understanding of the fundamental mechanisms governing 
the properties of the transition metals, the recently ob- 
served negative thermal expansion coefficient along the 
c-axis of one of these elements, the hep (a) phase of Ti [5j , 
stands out as an enigma. Especially since no other tran- 
sition metal so far has been shown to display such a be- 
haviour. 

The problem of finding connections between the elec- 
tronic structure of metals and alloys and peculiarities of 
their lattice properties has a long history, starting with 
the "third Hume-Rothery rule" concerning boundaries of 
phase stability in noble-metal alloys Q and its explana- 
tion by Jones in terms of touching of the Brillouin zone 
faces by the Fermi sphere [7] (for a review of further 
developments of these ideas, see Ref. [1]). The general 
concept of electronic topological transitions (ETT), in- 
troduced by L Lifshitz [9], that is, a coincidence of the 
Fermi level with a Van Hove singularity of the electronic 
density of states (DOS), is of crucial importance for un- 
derstanding these interrelations. Phase transitions and 
pre-martensitic anomalies of elastic moduli in alkali and 
alkaline-earth metals under pressure provide a clear ex- 
ample of the effects of the Van Hove singularities on the 
lattice properties 111- It turns out that the singu- 
larity in the electron DOS at the Fermi energy, N{Ep), 
should be visible also in elastic moduli and Debye tem- 
perature and, thus, in the thermodynamic properties 
of metals at low enough temperatures (the anomalies 
in phonon spectra with large enough wave vectors and 
thus in high-temperature thermodynamic properties are 
in general weaker, see Ref. Q and references therein). 



Since the thermal expansion is connected with the pres- 
sure derivatives of the elastic moduli, anomalies in the 
thermal expansion might be especially strong. In par- 
ticular, it can be proven thermodynamically that ETT 
in non-cubic metals should lead to a singular anisotropic 
thermal expansion at low enough temperatures 
The latter means that in principle it is always possible 
to prepare a textured material with zero thermal expan- 
sion. This conclusion, being interesting in itself, opens 
new ways to find nonmagnetic Invar systems. However, 
based on these general considerations alone it is impos- 
sible to predict the temperature region where the effect 
should be observable, or how far from the point of ETT 
the effect is still visible. Here we answer these questions 
based on direct microscopic calculations, in a framework 
of the density functional theory, and we address the re- 
cently discovered negative thermal expansion of a-Ti [sl . 

The occurrence of negative thermal expansion at low 
temperatures for non-cubic elemental solids have been 
known for quite some time, but only for elements outside 
the transition metal series. For instance , th e hexagonal 
close packed metals zinc and cadmium [13j have nega- 
tive thermal expansion coefficients along the a-axis {a±) 
for temperatures below ~ 75 K, while amongst the HIE 
group of the Periodic Table it is tin and indium that 
have negative thermal expansion coefficients along the 
basal plane (ax) and orthogonal to the basal plane (q;||), 
respectively [13]. Amongst the transition metals, Ti how- 
ever stands out. 

The analysis presented here is based on first principles 
density functional theory of the electron and phonon con- 
tributions to the total energy. We write the Helmholtz 
free energy as 
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where are the elastic constants, e the elastic strain, 
V the volume, F^^""^ the phonon free energy and F^^ the 
energy of thermal excitations in the electron subsystem. 
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In this expression the reference (zero) level is for a crys- 
tal at equilibrium conditions at zero temperature. The 



elastic constants were calculated from first principles [14 1. 

To evaluate the free energy contribution i^p'*°"^ which 
can be expressed as IgI 17| 



J dujg{uj,e)[— + kBTln{l 



)], (2) 



the phonon DOS g{Lo,e) has to be calculated. This was 
done within the quasi-harmonic approximation 



m 1 



where all anharmonic effects except the thermal expan- 
sion are neglected when calculating the temperature de- 
pendence of the phonons. In practice the phonon DOS 
was calculated by making small displacements of the 
atoms in a supercell (SC) 0. The directions of the 
displacements were [110] and [001] with amplitudes that 
were equal to ^0.4% of the lattice constant. The su- 
percell used was a 3x3x2 cell. Further details are found 
in RefllJ. The energy of thermal excitations of electron 
states i^^' was calculated by the standard expression by 
Sommerfeld and Frank [20]. 
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(3) 



where D{ep,e) is the calculated electronic density of 
states at the Fermi level. 

In Fig.l we compare our calculated phonon spectrum 
(at T = 0) with experimental values (at room tempera- 
ture). It is worthwhile to mention that a tight-binding 
calculation of the phonon dispersion for hep Ti has been 
published recently [21j, where parameters of the model 
were fitted to experimental data as well as to first prin- 
ciples calculations. The theoretical phonon dispersion 
curve in Fig. 1 agrees very well with the theoretical curves 



in Ref. [2l|. When comparing the theoretical and ex- 



perimental 2^ curves, we note an overall agreement, 
although certain differences can be identified. For in- 
stance, along the T — A direction the theory underesti- 
mates the frequencies in the lowest experimental branch, 
whereas the higher branches are reproduced with better 
accuracy, especially using the general gradient approxi- 
mation (GGA). Also, the calculated lowest branch along 
the T — K direction comes out somewhat too low com- 
pared to observations. The phonon DOS was then calcu- 
lated with the method of Ref. [2^ . 

By differentiating the free energy ^ with respect to 
£.„ = d{ln{V)) and Ec = d{ln{c/a)), it is possible to obtain 
an expression for the change in volume and structural 
property as a function of temperature. These changes 
are expressed in terms of equilibrium strains ,6^ and e[? 
at which dF/de^ = dF/dec = 0, and can be written in 
terms of the elastic constants and strain derivatives of 
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FIG. 1: The phonon dispersion of hep Ti at room tempera- 
ture and ambient pressure. The solid curve is the calculated 
frequencies from GGA and the dashed curve represents cal- 
culated frequencies from LDA. Both calculations are done at 
the experimental volume Vo — 15.9lA^. The filled circles are 
the experimental data of Stassis et al [2^ . 



the free energy 
where 
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Furthermore by differentiating ^ and ([5]) with respect 
to the temperature the following relations are obtained 
for the thermal expansion coefficients [12] 
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FIG. 2: Linear thermal expansion for hep Ti at ambient pres- 
sure. The solid lines are the theoretical calculation using GGA 
and the dashed line is from an LDA calculation. The filled 
circles are the experimental data of Nizhankovskii et al 
The filled squares are the experimental volume expansion co- 
efficient of Mal'ko et al [13]. The filled and empty triangle 
are the experimental data of aa and ac respectively of Pawar 
et al f25l|. The empty circles are the experimental data of 
Nizhankovskii et al scaled to give a volume expansion coeffi- 
cient in agreement with Refl24l. 



where a, = i^, a, = and 

By fitting free energies calculated at different strains 
and at a given temperature to polynomials of first degree 
in and second degree in Ec, the equilibrium strains 
can be obtained from Eqns. (|4]) and (O, and the thermal 
expansion coefficients can be calculated from Eqns. (fTO|l - 

In Figl2]we show the calculated thermal expansion co- 
efficients of a-titanium. The most important informa- 
tion to be extracted from this figure is that the observed 
negative thermal expansion coefficient along the c-axis 
is reproduced by our theory, where especially the calcu- 
lation based on GGA reproduce observations with the 
highest accuracy. It should be noted that GGA often is 
found to describe chemical bonding with better accuracy 
than LDA. The temperature interval for which ac is neg- 
ative is roughly 0-170 K, both in the observations and 
from the theory. The order of magnitude of ac and aa 
is also the same when comparing experiment and theory. 
Figure also shows that theory reproduces, with good 
accuracy, the volume expansion coefficients of Ref. |24, 
especially at somewhat elevated temperatures. We also 
note that based on thermodynamic relations /3 should 
approach zero at T = K, which our theoretical curves 
do. 



The fact that both the measured and calculated ther- 
mal expansion coefficients along the c-axis of Ti are nega- 



tive at low temperatures strongly suggests the uniqueness 
of elemental Ti among transition metals, although the ab- 
solute value of the measures low temperature expansion 
coefficient is still somewhat uncertain. The measured 
data of Ref. |5| (filled circles) have in Fig. [5] been scaled 
(open circles) to reproduce room temperature values of 
(3, aa and ac, and it is found that these scaled values 
compare better with our theory (Fig. Although a 
slight calibration error in Ref. i can not be excluded, 
there is good reason to view the negative value of ac at 
low temperatures as a true materials property of a-Ti. 

As we will show below the microscopic origin of the 
negative thermal expansion for ac of Ti is due to the 
closeness to a saddle point van Hoove singularity of the 
electronic structure. To illustrate this singularity we pro- 
ceed with an analysis of the Fermi surface. In order to 
do this we show in Fig. [3] the calculated Fermi surface at 
the equilibrium volume for three different values of the 
out-of-plane lattice constant, c. The figure shows that as 
the c lattice constant decreases the inner ellipsoidal sur- 
face at the F-point and the Fermi surface centered at the 
A-point, become connected along the TA line. The elec- 
tronic structure as revealed by the Fermi surface shown 
in Fig. 3 thus demonstrates the presence of a saddle point 
Van Hove singularity, which is associated with a singular 
contribution to the density of states at the Fermi level 

N{Ef), 6N{Ef) VEf - Eci9{Ef - Ed) , where 

9{x) fulfills: 9{x > 0) = ^9{x < 0) = and Ed is the 
critical point energy The energy difference be- 

tween Ef and the energy of the critical saddle point, at 
the theoretical equilibrium volume and a c=co, has been 
calculated to be Ef — Ed ~ 44 meV. Another critical 
point, associated with the appearance of a new ellipsoid 
around the K symmetry point (not shown in Fig. [3] ) 
has been found in the calculations, giving rise to the sin- 
gular contribution 5N{Ef) VEc2 - Ef9{Ec2 ~ Ef) 
to N{Ef)- However since Ec2 — Ef 142 meV and 
\d{Ec2-EF)/de\ < \d{EF - Ed)/ de\, it is clear that the 
saddle-point topological transition, at Ed, gives rise to 
the strongest singular contribution to N{Ef)- 

By calculating the derivatives of N{Ef) with respect 
to the two different types of strains we have found that 

dN{EF)/dey - 0.75 eV'^ and dN{EF)/dec 0.77 

eV~^. Since the singularities in N{Ef) influence the elas- 
tic moduh, thus effecting the Debye temperature Q, it is 
natural to attribute the main reason for the anisotropic 
thermal expansion in a— titanium to the saddle-point 
Van Hove singularity near the Fermi level. 
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FIG. 3: Calculated Fermi surface of hep Ti at T=0 equi- 
librium volume for three different lattice constants c, at 
c = 0.988co (a), at c = co (b) and at c = 1.012co (c). Here co 
corresponds to the (T=0) equilibrium lattice constant. 
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